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We reexamine the two-point correlation of density maxima in Gaussian initial conditions. Spatial 
derivatives of the linear density correlation, which were ignored in the calculation of Bardeen et 
al. [Astrophys. J. , 304, 15 (1986)], are included in our analysis. These functions exhibit large 
oscillations around the sound horizon scale for generic CDM power spectra. We derive the exact 
leading-order expression for the correlation of density peaks and demonstrate the contribution of 
those spatial derivatives. In particular, we show that these functions can modify significantly the 
baryon acoustic signature of density maxima relative to that of the linear density field. The effect 
depends upon the exact value of the peak height, the filter shape and size, and the small-scale 
behaviour of the transfer function. In the ACDM cosmology, for maxima identified in the density field 
smoothed at mass scale M « 10 12 — 10 14 Mg/li and with linear threshold height v — 1.673/cr(M), 
the contrast of the BAO can be a few tens of percent larger than in the linear matter correlation. 
Overall, the BAO is amplified for v J>1 and damped for v <Cl. Density maxima thus behave quite 
differently than linearly biased tracers of the density field, whose acoustic signature is a simple scaled 
version of the linear baryon acoustic oscillation. We also calculate the mean streaming of peak pairs 
in the quasi-linear regime. We show that the leading-order 2-point correlation and pairwise velocity 
of density peaks are consistent with a nonlinear, local biasing relation involving gradients of the 
density field. Biasing will be an important issue in ascertaining how much of the enhancement of 
the BAO in the primeval correlation of density maxima propagates into the late-time clustering of 
galaxies. 
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I. INTRODUCTION 

Sound waves propagating in the primordial photon- 
baryon fluid imprint a oscillatory pattern in the 
anisotropies of the Cosmic Microwave Background 
(CMB) and in the matter distribution, whose charac- 
teristic length scale r s is the sound horizon at the re- 
combination epoch r s « 105 /i _1 Mpc for the cur- 
rently favoured cosmological models. While experiments 
have accurately measured this fundamental scale and 
its harmonic series in the temperature and polarisation 
power spectra of the CMB, this acoustic signature has re- 
cently been detected in the correlation function of galax- 
ies 9- There is also weak evidence for the baryon 
oscillations in the correlation function of clusters [H . In 
the 2-point correlation, the series of maxima and minima 
present in the power spectrum translates into a broad 
peak at the sound horizon scale. Since the latter can 
be accurately calibrated with CMB measurements, the 
baryon acoustic oscillations (BAO) have emerged as a 
very promising standard ruler for determining the angu- 
lar diameter distance and Hubble parameter [5J . Measur- 
ing the BAO at different redshifts thus offers a potentially 
robust probe of the dark energy equation of state. 

In linear theory, the amplitude of the baryon acoustic 
peak increases while its shape and contrast remain un- 
changed. However, the clustering of galaxies does not 
fully represent the primeval correlation. Mode-coupling, 
pairwise velocities and galaxy bias are expected to al- 
ter the position and shape of the acoustic peak and, 
therefore, bias the measurement @. The evolution of 



the acoustic pattern in the 2-point statistics of the mat- 
ter, halo or galaxy distributions has been studied us- 
ing both numerical simulations @, Q and analytic tech- 
niques based on the halo model or perturbation the- 
ory [T3, E3, E3, H^> E3, HM> H^> E3, H3> IS H3- 

Yet the results of these studies do not always agree and 
the impact of nonlinearities on the matter and galaxy 
power spectrum remains debatable. For instance, Ref- 
erences [ill, E3 argue that any systematic shift (i.e. 
not related to random motions or biasing) must be less 
than the percent level owing to the particularly smooth 
power added by nonlinearities on those scale, and to 
the cancellation of the mean streaming of (linearly) bi- 
ased tracers at first order. On the other hand, Refer- 
ences [HI, Oj], [H, G3 have shown that mode-coupling 
modifies the acoustic pattern in the correlation of dark 
matter and haloes, and generates a percent shift towards 
smaller scales. Despite their redshift dependence, these 
shifts appear to be predictable and could be removed 
from the data 20 . 

There is a broad consensus regarding the shape of the 
acoustic peak. In light of the nonlinear gravitational evo- 
lution of matter fluctuations, it is sensible to expect a 
baryon acoustic peak less pronounced in the late-time 
clustering of galaxies than in the linear theory correla- 
tion. This can be shown to hold for any local transfor- 
mation of the density field [H, [H, [24[. Such biasing 
mechanisms do indeed predict a damping of the baryon 
acoustic features in the 2-point statistics of the galaxy 
distribution [1, [ij Galaxies, of course, form a dis- 
crete set of points but one commonly assumes them to 
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be a Poisson sample of some continuous field. Still, the 
extent to which those models are an accurate approxi- 
mation to the clustering of galaxies remains unclear. No- 
tice also that a reconstruction of the primordial density 
field could significantly restore the original contrast of 
the acoustic oscillation [221 ] . 

The main objective of this paper is to demonstrate that 
the BAO in the correlation of tracers of the density field 
can be noticeably modified if we consider local biasing 
relations more sophisticated than local transformations 
of the density field [25j, |26j. To this purpose, we will 
examine the clustering of density maxima in the initial 
cosmological density field. In this respect, we will assume 
that the initial fluctuations are described by Gaussian 
statistics. This assumption is remarkably well supported 
by measurements of the CMB and large-scale structures 
[23,[2^]. Density peaks form a well-behaved point-process 
whose statistical properties depend not only on the un- 
derlying density field, but also on its first and second 
derivatives. Therefore, although the number of density 
maxima is modulated by large-scale fluctuations in the 
background, their clustering properties cannot be derived 
from a continuous field approach in which the peak over- 
density would depend upon the value of the matter den- 
sity only. Interestingly however, we shall see that, at 
large separations, the peak correlation and pairwise ve- 
locity are consistent with a nonlinear biasing relation in- 
volving gradients of the density field. 

In a seminal paper, Bardeen, Bond, Kaiser & Sza- 
lay (hereafter BBKS) [2!| provided a compact ex- 
pression for the average number density of peaks in a 
three-dimensional Gaussian random field, etc. Further- 
more, they obtained a large-scale approximation for the 
correlation function of peaks which, at large threshold 
height, tends toward the correlation of overdense re- 
gions [3(1 HHHl] as ^ should be. However, BBKS deter- 
mined the correlation function of density maxima only in 
the limit where derivatives of the 2-point function of the 
density field can be ignored. As we will see below, these 
correlations can greatly influence the large-scale correla- 
tion of density maxima for generic Cold Dark Matter 
(CDM) power spectra. It is also worth noticing that 
the statistics of Gaussian random fields in a cosmolog- 
ical context has received some attention in the litera- 
ture [H, [H, Hil [H, H3] • Some of these results have been 
applied to the mass function and correlation of rich clus- 
ters for example [H, [H, E3 ■ The present work mainly 
follows the analytic study of BBKS, and the lines dis- 
cussed in [4l|, [43] , where 2-point statistics of the linear 
tidal shear are investigated. We refer the reader to [43| 
for a rigorous introduction to the statistics of maxima of 
Gaussian random fields. 

The paper is organised as follows. Section |TT] intro- 
duces a number of useful variables and correlation func- 
tions. Section [TTT1 is devoted to the derivation of the lead- 
ing order expression for the large-scale asymptotics of the 
peak correlation. Our result can be thought as arising 
from a specific type of nonlinear local biasing relation 



including second spatial derivatives of the density field. 
In Sec. HV1 we explore the impact of these derivatives 
on the amplitude and shape of the correlation of density 
maxima. Our attention focuses on the baryon oscillation, 
across which the amplitude of the linear matter correla- 
tion varies abruptly. It is shown that the BAO of density 
maxima can be amplified relative to that of the mat- 
ter distribution. Section [V] deals with the peak pairwise 
velocity. Its leading order contribution is found to be 
consistent with the nonlinear local bias relation inferred 
from the 2-point correlation of peaks. A final section 
summarises our results. 



II. PROPERTIES OF COSMOLOGICAL 
GAUSSIAN DENSITY FIELDS 

We review some general properties of Gaussian random 
fields and provide explicit expressions for the correlations 
of the density and its lowest derivatives. We show that 
the latter are not always negligible in CDM cosmologies. 



A. Useful definitions 

We will assume a ACDM cosmology with normalisation 
amplitude <7g — 0.82, and spectral index n s = 0.96 [HI]. 
The matter transfer function is computed using publicly 
available Boltzmann codes [3] . The position of the BAO 
in the linear matter correlation function is close to ss 
105.0 /i^Mpc. 

Let q designate the Lagrangian coordinate. We are 
interested in the three-dimensional density field S(q) and 
its first and second derivatives. It is more convenient to 
work with the normalised variables v = S(q)/ao, r\i = 
di8{q)/(j\ and (ij — didj5(q)/a 2 , where the <jj are the 
spectral moments of the matter power spectrum, 

/>oo 

a) = J dln/fc/fc 2j ' A 2 (k) . (1) 

A 2 (k) = A 2 (k)\W(k, Rf)\ 2 denotes the dimensionless 
power spectrum of the density field smoothed on scale 
Rf with a spherically symmetric window W(k, Rf). 

The best choice of smoothing is open to debate. 
Among the popular window functions, the top hat filter 
is compactly supported and has a straightforward inter- 
pretation within the spherical collapse model. Notwith- 
standing this, oscillations that arise in Fourier space do 
not lead to well defined spectral moments a.j with j >2 
for CDM power spectra. This can be understood by ex- 
amining the high-fc tail of the CDM transfer function. 
Neglecting the baryon thermal pressure on scale less than 
the Jeans length, the small-scale matter transfer func- 
tion behaves as T(k) cx ln(1.8fc)/fc 2 [Hi EH, which clearly 
leads to divergences when the integer j is larger than one. 
By contrast, a Gaussian window function ensures the 
convergence of all the spectral moments for any realistic 
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FIG. 1: A comparison between the cross-correlation of the 
density field, £(r), and that of its first and second derivatives, 
E(r) and ijj(r) respectively (see eq. [6]). Results are shown 
as a function of the Lagrangian separation r for the ACDM 
cosmology considered in the present work. The density field 
is smoothed with a Gaussian filter of characteristic scale 
Rf = 5 /i _1 Mpc (i.e. a mass scale Mf = 1.5 x 10 14 M Q /h). 
Dashed lines denote negative values. All the correlations are 
normalised to unity at zero lag. 



matter power spectra. Consequently, we shall mostly rely 
on the Gaussian filter throughout this paper, although 
the top hat filter will also be considered briefly in Sec. 
IIVI Note that a Gaussian filter of characteristic width 
Rf encloses a mass Mf = (27r) 3 / 2 pRj a few times larger 
than that encompassed by a top hat filter of identical 
smoothing radius. 

Following BBKS, we also introduce the parameters 7 = 
<t 2 /(<7o<t 2 ) and R± = \fZo\lu-x for subsequent use. The 
spectral width 7 reflects the range over which A 2 (fc) is 
large, while i?* characterises the radius of peaks. For the 
special case of a powerlaw power spectrum with Gaussian 
filtering on scale R / , these parameters are given by = 
(n + 3)/(n + 5) and i? 2 = 6R 2 f /(n + 5). For CDM power 
spectra, 7 ~ 0.5 — 0.7 when the smoothing length varies 
over the range 0.1 <Rf <10 /i _1 Mpc. 



B. Correlation of the density and its derivatives 

Calculating the 2-point correlation of density peaks re- 
quires knowledge of the auto- and cross-correlations of 
the various fields. These objects can be decomposed into 
components with definite transformation properties un- 
der rotations. Statistical isotropy and symmetry implies 
that, in position space, the most general ansatz for the 
isotropic sector of the 2-point correlations of these fields 



reads 

M9iM«fa)> = £M (2) 

M<?iM<72)) = 3(r)fi 
(Vi(li)Vj(92)) = Si(r)nfj + £ 2 (r)(5.y 

= Ili(r) nflfm + n 2 (r) (hShn + rtfim + fmkl) 
(Cy(<Zl)Czm(9 2 )) = ^l(r)fifjflf m 

+ * 3 (r) (r + rif m 5ji + fjfid 

im T" ^j^m^il 

) + *s(r) (SijSim + 6aSjm + km^n) , 

where r = \q 2 — q±\ is the Lagrangian separation, f» = 
Ti/r and the functions £, S, Ej, IT and 'F; depend on 
r only. We emphasise that these correlation functions 
transform as scalar under rotations. Note also that these 
expressions are valid for any arbitrary random field. For a 
cosmological Gaussian density field however, these func- 
tions can be summarised as follows : 
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je(x) arc spherical Bessel functions of the first kind. In 
the limit r —> 0, all the correlation functions vanish but £, 
E 2 and ^5, which tend towards 1, 1/3 and 1/15, respec- 
tively. Averaging over the direction r of the separation 
vector thus yields 



S(r) 
— Sii 



■^JdCl i .(r] i (q 1 )r] j (q 2 )) = 

1 /" ip(r) 
— J dfif. (Cij(qiKim{q 2 )) = -j^- 



(4) 
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for the covariances of the fields T)i and Qj , where we have 
defined 

E(r) = S!(r)+3E 2 (r) (5) 
i>{r) = *i(r) + 10* 3 (r)+15* B (r). 

The angular average of the other correlation functions 
vanishes, except that of the density correlation of course. 

E(r) and tf>(r) can be expressed in terms of the 
derivatives of the density correlation using relations like 
(ViVj) — ~didj^(r) etc. For a density correlation that 
falls off as a powerlaw r~ n ~ 3 , E(r) and ip(r) decay 
as r~™~ 5 and r~™~ 7 , respectively. This derivation as- 
sumes a powerlaw power spectrum with a fair amount 
of power at short wavenumbers. Hence, as recognised in 
BBKS, neglecting the derivatives of the density correla- 
tion should be a reasonable approximation when n < — 1. 

This simple argument may not hold for CDM cosmolo- 
gies since the index n is a smooth function of the separa- 
tion r. Namely, it is n ~ —2 when r ~ 10 /i _1 Mpc, and 
increases to attain a value of the order of unity on scale 
r >60 h M.pc. For illustration purpose, the functions 
£, E and i/j are plotted in Figures □ and [U for the ACDM 
cosmology considered here. The filtering length is Rf = 5 
and 1 /j _1 Mpc, respectively (The reason for choosing 
these values will become apparent below). Retaining 
only the density correlation appears to be a good ap- 
proximation on scales larger than a few smoothing radii. 
However, the relative amplitude of the cross-correlations 
strongly depends upon the filtering scale. Namely, both 
E(r) and ip( r ) increase relative to £(r) with increasing 
smoothing length. Yet another striking feature of Fig- 
ures [T] and [2] is the oscillatory behaviour of E(r) and 
ij){r). The large oscillations are caused by rapid changes 
in the linear matter correlation across the baryon acous- 
tic peak. Notice that both E(r) and tp(r) are positive 
at distances r w 100 — 110 /i _1 Mpc. On these scales, 
when Rf = 1 /i _1 Mpc, E(r) reaches to 3 per cent of the 
density correlation while ip(r) is negligible. At the large 
smoothing length however, they nearly reach 20 and 10 
per cent of the density correlation, respectively. 

These results suggest that, for generic CDM power 
spectra, the derivatives of the density correlation could 
have a significant impact on the correlation of density 
maxima, especially in the vicinity of the baryon acoustic 
feature. This motivates the calculation presented in the 
next Section. 



III. CORRELATION OF DENSITY MAXIMA 

Owing to the constraints on the derivatives of the den- 
sity field, calculating the n-point correlation function of 
peaks requires performing integration over a joint prob- 
ability distribution in lOn variables. Therefore, even the 
evaluation of the 2-point correlation of density maxima 
?pk(0 proves difficult. Here, we derive the leading order 
expression that includes, in addition to the linear matter 
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FIG. 2: Same as Fig. [T] but for a smoothing length Rf = 
1 h~ 1 Mpc (M f = 1.2 x 10 12 M Q //i). The correlation function 
ip(r) (not shown) is less than 10 -6 at distances larger than 
^30 /i _1 Mpc. 



correlation £(r), the contribution of the angular average 
functions E(r) and ip(r). We also show that the large- 
scale asymptotics of the peak correlation can be thought 
as arising from a specific type of nonlinear biasing rela- 
tion involving second derivatives of the density field. 



A. The Kac-Rice formula 

As shown in BBKS for instance, the correlation of den- 
sity extrema (maxima, minima and saddle points) can 
be entirely expressed in terms of 6(q) and its derivatives, 
f]i(q) and Cijio)- I n the neighbourhood of an extremum, 
the first derivative r\i is approximately 

r)i(q) w V3R- 1 tij(Qp) («~ %) ■ ( 6 ) 

3 

Using the properties of the Dirac delta, the number den- 
sity of extrema can be written as 

n ext (g) = J2 S3 (l- Q P ) = ^ |detC(g)|<5 3 [77(g)] , (7) 

provided that the Hessian is invertible. The delta 
function <5 3 [77] ensures that all the extrema are included. 
In this paper however, we are interested in counting the 
maxima solely. Consequently, we further have to require 
djilp) be negative definite at the extremum position q p . 
Note that, later, we will also restrict the set to those 
maxima with a certain threshold height. The average 
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number density of maxima eventually reads 

(nMQ)) = ^-(\detaq)\S 3 [v(q)])- (8) 



This expression, known as the Kac-Rice formula 43, 46, 
S3, 53, S§| , holds for arbitrary smooth random fields. In 
the general case of a random field in N dimensions, it is 
trivial to show that the mean density of maxima scales 

as (« P k) fx R± N « RJ N ■ 

The 2-point correlation function of density peak is de- 
fined such that 



i + £ P k(» = (n P k(gi)"pk(g 2 ))/(«pk> 



(9) 



is the joint probability that a density maxima is in a 
volume dVi about each Let X be the diagonal matrix 
of entries diag(a;i, x%, x%) where x\ > x 2 > x 3 is the 
non-increasing sequence of eigenvalues of the symmetric 
matrix — £. The condition that the extrema are maxima 
implies x 3 > 0. Therefore, the correlation function of 
peaks is given by 



(10) 



M 2 Rl 

3 3 



(|detCi||det&| e(x 3 )8(y 3 ) 



di/Ki 6 Cidi/ 2 d 6 C 2 |detCi||detC 2 | 0(x 3 )0(y 3 



M 2 R$ 

x P(r)x = 0,i/i, (i,m = 0,v 2 ,b;r) , 

where, for shorthand convenience, subscripts denote 
quantities evaluated at different Lagrangian positions, 
and d 6 £ = Yli<j ^-dj is the usual Lebesgue measure 
on the six-dimensional space of symmetric matrices. 
Here and henceforth 6(x) designates the Heaviside step- 
function, i.e. 6{x) = 1 for x > and zero otherwise. 



B. Two-point probability distribution 

The joint probability distribution of the density 
fields together with its first and second derivatives, 
P (r)i , v\ , Ci , r)2 , vi , £ 2 ; r) , is given by a multivariate Gaus- 
sian whose covariance matrix C has 20 dimensions. This 
20 x 20 matrix may be partitioned into four 10 x 10 block 
matrices, M = (y^J) = (y 2 y 2 ) m the top left corner 
and bottom right corners, B = (j/iJ/ 2 ) and its trans- 
pose in the bottom left and top right corners, respec- 
tively. The components (a, A — 1, ...,6 of the ten- 
dimensional vector y T = (r^, v, £4) symbolise the entries 
ij = 11,22,33,12,13,23 of Cij- To emphasise that the 
entries (a transform as a tensor under rotation, we shall 
also label them as the matrix £ in what follows. 

The matrices M and B can be further decomposed into 
block sub- matrices of size 4 and 6, 



M=| Ml M 3 
1 M 3 M 2 



B = 



Bi BJ 
B3 B 2 



(11) 



Unlike the M, which describe the covariances at a single 
position, the matrices B^ generally are functions of the 
separation vector r. Using the harmonic decomposition 
of the tensor products • • • ® f, they can be written as 



Bj(r) = Br+^Bf>)>™ 

1=1 

Mr) = B«'°+^B^(r)y™(f) 
l=\ 

B a (r) = B° 3 '°+J2B e r(r)Yr(r). (12) 

Yg m (f) are spherical harmonics and B^' m (r) are matrices 
which satisfy (B^ m )t = (-l) m Bf m . Only multipoles up 
to £ = 4 appear in the harmonic decomposition since the 
correlations given in eq. ([3]) involve products of up to four 
unit vectors f. The monopole terms are 
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0,0 _ ( E(r)/3I 3 xi 
Oixs t{r) 

R o,o _ ( VW/15A 3x3 

2 " ^ 3x3 ^(r)/15I 

R o,o _ ( 3x3 -7S(r)/31 3x i 

3 I 3x3 3x i 



(13) 



where 




(14) 



I is the 3x3 identity matrix and li X 3 = (1,1,1) etc. 
The matrices Mj are readily obtained as M, = B^' (0). 
An explicit computation of the higher multipole matrices 
is unnecessary here as we confine the calculation to the 
monopole contribution. 

It is important to note that the joint density 
P(y 11 y 2 ,r) preserves its functional form under the action 
of the rotation group SO (3). However, in a given frame of 
reference, P(yi, y 2 , r) does change when f moves on the 
unit sphere. Does this mean that £ p k(»") truly depends 
on the direction of the separation vector r ? No, as it 
should be clear from eq. (fTTjl where the volume measure 
|det£|d 6 C is a rotational invariant. More precisely, the 
volume element d 6 £ can be cast into the form 



d 6 C = 8tt 2 |A(a;)|d 3 a;dR 



(15) 



where the XiS are, as before, the three ordered eigenvalues 
of — £, d 3 ir = dxidcc 2 da; 3 and A(x) = Ilfo-OEi — x j) is 
the Vandermonde determinant. dR is the Haar measure 
(for the Euler angles for example) on the group SO (3) 
normalised to J dR = 1 . The peak correlation thus is 
proportional to 



J dRidR 2 P(?7i = 0, ui, Ci,r/a = 0, u 2 , (2, r) 



(16) 
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where the integral runs over the two SO (3) manifolds that 
define the orientation of the principal frames of ^ and £2 
relative to the frame of reference. Alternatively, we can 
choose the coordinate system such that the coordinate 
axes are aligned with the principal axes of £1. In this 
new coordinate frame, the above integral becomes 



(17) 



where R is an orthogonal matrix that defines the orien- 
tation of the eigenvectors of £2 relative to those of £1. 
This demonstrates that only the monopolc component of 
P(Ui, y 2 , r ) contributes to the peak correlation function. 
Therefore, ^pk(^) is invariant under rotations of the coor- 
dinate system, namely, it is a function of the separation 
r only. 



tribution M, e.g. |B| < M. Following j4l|, [42j], the 
quadratic form which appears in the probability distri- 
bution P(y 1 ,y 2 ;r), 



P{ynV2\ r ) 



(27r) iu IdetCI 1 ^ 



(19) 



where detC « |detM| 2 = 4 2 (1 - 7 2 ) 2 /(15 10 3 8 ) is the de- 
terminant of the covariance matrix C, can be computed 
easily using Schur's identities. Expanding the exponen- 
tial in the small perturbation B yields, to first order, 



e -Q(yi-y 2 ,r) ~ (1 + yjM^BM^y^ e^C^te) , (20) 



C. Large scale asymptotics 



where the quadratic form Q(y 1 , y 2 ) can be recast as 



To obtain the correlation function of peak, we need first 
to calculate the 2-point probability distribution function 
averaged over the unit sphere for the variables y T = 

P(y 1 ,y 2 ,r) = ^fdn t .P(y 1 ,y 2 ,r). (18) 

In the large-distance limit (r 3> 1), the cross-correlation 
matrix is small when compared to the zero-point con- 



2Q = v^ (™ +tr &) 



1 



3tr(C 1 2 )-(trCi) 2 J+1^2, 
(21) 

in agreement with the results of BBKS. The calculation 
of yjM~ 1 BM~ 1 y 2 is tedious but straightforward. Fortu- 
nately, only the monopole terms B°'° survive after aver- 
aging over the directions f. After further simplification, 
the result can be reduced to the following compact ex- 
pression : 



■ fdn~ r yjM- 1 BM- 1 y 2 = | [3tr (C1C2) - tr& tr£ 2 ] *l>(r) + {tr£x tr£ 2 ty(r) + 7 2 £(r)] (22) 

+ [£(r) + - 2 7 2 (trCi trC 2 + v x v 2 ) E(r) + 7 + mtrCi) [£(r) + i>(r) - (l + 7 2 ) E(r)] } (l - 7 2 )~ 2 . 

I 



The invariance under rotation requires that P(yi, 2/2^) 
be a symmetric function of the eigenvalues, and thus a 
function of tr (Ci C2) > M = 0, 1, . . . 

Since the above expression depends only upon the rel- 
ative orientation of the two principal axes frames of £1 
and £2 (through the presence of tr(£i£ 2 )), we choose a 
coordinate system whose axes are aligned with the prin- 
cipal frame of £1. With this choice of coordinate, we 
define £1 = — X and £2 = — RYR T , where R is an or- 
thogonal matrix that defines the relative orientation of 
the eigenvectors of £2. X and Y are the diagonal ma- 
trices consisting of the three ordered eigenvalues Xj and 
Hi of the Hessian —didjV. The properties of the trace 
imply that tr£j = — trX, tr(£ 2 ) = tr(X 2 ) (and similarly 
for £ 2 ), while the term tr(£i£ 2 ) = tr(XRYR T ) depends 
explicitely on the rotation matrix R. 



The integral over the SO (3) manifold that describes the 
orientation of the orthonormal triad of £1 is immediate. 
The result is 27r 2 (and not 87r 2 ) as we don't care whether 
the axes are directed towards the positive or negative 
direction. The integral over the second SO (3) manifold 
involves 

/ dRtr (XRYR T ) = \ trXtrY , (23) 
iso(3) 3 

and yields cancellation of the first term in the right- 
hand side of eq. (fTT)) . To integrate over the eigenvalues 
of £1 and £2, we transform to the new set of variables 
{ui, Vi, w iy i = 1, 2}, where 

U\ = X\ + X 2 + X 3 

Vi = (sci - x 3 ) /2 

id = (xi - 2x2 + ar 3 ) /2 . (24) 
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The variables (1/2, 1)2,1112) are similarly defined in terms 
of the yi. We will henceforth refer to u as the peak cur- 
vature. 

Our choice of ordering imposes the constraints Vi > 
and — vi < Wi < «». The condition that the density ex- 
trema be maxima, i.e. all three eigenvalues of the Hessian 
Qij are negative, translates into [iii +Wi) > 3^. Another 
condition, m > 0, should also be applied if one is inter- 
ested in selecting maxima with positive threshold height. 

For shorthand convenience, and to facilitate the com- 
parison with the calculation of BBKS, we introduce the 
auxiliary function 



F(ui,vi,wi) 



|detX| A(x) 



(25) 



(iti — 2wi) {u\ + wi) 2 — 9vf v\ (y\ — w 2 ) 



F(ui, vi,w\) measures the degree of asphericity expected 
for a peak and can be used to determine the proba- 
bility distribution of ellipticity v\/u\ and prolateness 
wi/ui . It scales as oc in the limit Ui 3> 1. 



D. The peak correlation £ pk (^, r) 

For sake of generality, we will present results for the 
cross-correlation ^ p k(i / i, ^2, f) between two populations 
of density maxima v\ ^ v 2 identified at smoothing scale 
R\ 7^ R2- However, we shall focus shortly on the auto- 
correlation $pk(^, t) [y\ = v-i = v and R\ = R 2 = Rf), 
which is more directly related to the clustering proper- 
ties of dark matter haloes of a given mass or galaxies and 
clusters of a given luminosity spanning a narrow redshift 
range. It may also be interesting to work out the cor- 
relation of peaks with a fixed height but identified at 
smoothing radii R> Rf, which can be thought as mim- 
icking the statistical properties of haloes above a given 
mass. However, we will not consider this correlation here 
since it requires a solution to the cloud-in-cloud prob- 
lem [Hoj at the location of density maxima. 

Let n pk = n p ^{i/) hereafter denote the differential den- 
sity of peaks in the range v to v + dv. The expectation 
value of the product of the local peak densities that ap- 
pears in eq. ([9|) is then 



£ P k {1/1,1/2, r) = 



1 5 b 3 



5q4 



R 



-6 



1 / II {duidv i dw l F{u i ,v l ,w i )} ^ {i/ 1: v 2 ,ui,u 2 , 



r) e 



(26) 



W 2 (2,) 6 * 
where 

$oOi, V2, ux,u 2 ,r) = {uiu 2 [ip{r) + 7 2 £M] (27) 

+ V1V2 [£(r)+ J 2 ip{r)] -2 1 2 {u 1 u 2 + v lV2 )Y I {r)- 1 {u 1 v 2 +u 2 vi) [f(r) + i/)(r) - (l+ 7 2 )S(r)]} (l- 7 2 )~ 2 , 

I 



is equation (|17p averaged over the relative orientation 
of the frames spanned by the eigenvectors of £1 and (2 ■ 
<&o depends on the separation r through the correlation 
functions £(r), S(r) and ijj{r) only. Furthermore, the 
quadratic form Q simply is 



2Q = ui 



{ui - JUl) 
1 - 7 2 



15vf + 5w 2 + 1 <-> 2 



(28) 



in the variables 

The integration over the variables Vi and tu, is lengthy 
but straightforward. We refer the reader to BBKS for 
the details since the calculation now proceeds along sim- 
ilar lines. Let us mention that the allowed domain of 
integration is the interior of a triangle bounded by the 



points (0,0), {ui/4, — Uj/4) and {ui/2, Ui/2). As shown 
in BBKS, the differential density of peak of height v can 
be cast into the form 



n p k{i/) = 



1 



WRl 



- 2 / 2 Go( 7 ,7^) 



(29) 



where Go is the zeroth moment of the peak curvature u. 
Higher moments are written in explicit compact form as 



Gn{l,u) 



dxx n f{x) 



^(1"7 2 ) 



(30) 



Using this result, the correlation of peaks can be rear- 
ranged as follows : 



£ P k(fi,^2,r) = G (7,7^i) ^0(7,7^2) 1 Yi S du */K 



, 2 , AFf) 
1 



$ {i/i,i/ 2 ,ui,u 2 ,r) . (31) 



For sake of completeness, 



1 



as demonstrated in BBKS, who noted also that the 













j^Erf 




+ Erf 




j 








V22 





(32) 
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asymptotic limits of this function include a cancellation 
to eighth order at small x, and the x 3 law expected for 
density maxima at large x. 

The integration over x must generally be done numer- 
ically. It is worth noticing that, while the exponential 
cxp[— (x — w) 2 /2(l — 7 2 )] decays rapidly to zero, x n f(x) 
are monotonically and rapidly rising. As a result, the 
functions G„(7, w) are sharply peaked around their max- 
imum. For large values of u>, we find that Go and G\ 
asymptote to 

G (7,w) « w 3 - 37 2 c J + S (7)^ 2 e- A(7) " 2 (33) 
Gi(7,w) m w 4 + 3w 2 (1 - 2 7 2 ) +B 1 ( 7 )w 3 e-^^ 2 . 



The coefficients A(j), -60(7) and -61(7) are obtained from 
the asymptotic expansion of the Error function that ap- 
pears in eq. ([33]). We have explicitly 



A = 



5/2 



5 7 2 ) 



j Bo = 



432 



r, B 1 = 



4S n 



lM9 - 5 7 2 ) 5/2 ' (9" 5 7 2 ) 

(34) 

The rest of the calculation is easily accomplished. The 
2-point correlation function of peaks eventually reads 



£ P k(^i,^2,7-) = {{v\ ~1 u i) ( v i -ju 2 ) £(r) 

+ (ui - 7^1) (u 2 - 7^2) ?/K r ) - [(^1 - 7"i) (7^2 - u 2 ) + (7^1 - ui) (u 2 - 7S2)] 7 s ( r )} (l - 7 2 ) 



(35) 



where we have introduced the mean curvature 72(7, 71/) — 
Gi/Gq. Also, the notation is such that Ui = ^(7,72/,-). 
The function 11(7,71/) is accurately fitted by eq. (4.4) 
of BBKS, which is constructed to match the asymptotic 
large v expansions of Go and G\ given in eq. (|34|) . In the 
special case V\ = v% = j/, the 2-point correlation of peaks 
simplifies to 

£ pk (i/, r) = 6 2 (1/, 7) £(r) + 6,(1/, 7) S(r) + 6 2 (i/, 7) ^(r) , 

(36) 

where the bias functions 6„, b^ and 6^ are 



b c (u, 7) 
M 1 ') 7) 



v 



I-7 2 
u — 71/ 



1 - 7 2 

2jb u (i/, 7) &c(j/,7) 



(37) 



The sign convention is chosen such that all three bias 
parameters are positive when v — > 00. Notice that b v is 
precisely the amplification factor found by BBKS when 
derivatives of the density correlation function are ne- 
glected. 

Equation (f3"6")l . which holds for any value of the peak 
height v and the smoothing length Rf , is the main result 
of this Section. It describes the asymptotic behaviour 
of the peak correlation function in the limit where the 
correlation functions £(r), S(r) and tp(r) are much less 
than unity. 




peak height v 



FIG. 3: Bias factors bl(v, 7), 677(1/, 7) and 6?(i/, 7) as a func- 
tion of the peak height v. The density field is smoothed on 
scale Rf = 5 /i -1 Mpc with a Gaussian filter. This leads to a 
correlation strength 7 = 0.676. Dashed curves indicate nega- 
tive values. The dotted curves are the asymptotic expansions 
given in eq. 



E. The bias parameters 

To gain some insight into the behaviour of the peak 
correlation function £pk(^, r), we have plotted in Fig. [3] 
the biasing parameters b v , b v and b^ as a function of 



the peak height. Again, the density field is smoothed on 
scale Rf — 5 h~ 1 Mpc with a Gaussian filter. The dotted 
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curves show the following large v approximations, 
K{v>l) ~ v- - 



b c {v,j) w — 



6,(1/, 7) « 6 1- 



(38) 



obtained from the asymptotic expansions of Go and G\ 
(eq. l34|) . They provide a good match to the bias param- 
eters when the peak height is larger than ~ 2. As we 
can see, b n tends towards the constant value of 6 when 
v — > oo. Moreover, for a threshold height less than unity, 
b v is negative and of absolute magnitude larger than b 2 ,. 
This is also true in the intermediate region v ~ 1 — 2. For 
these threshold heights, both b v and b n vanish while the 
bias parameter b^ is of the order of a few. Consequently, 
the correlation of density maxima, albeit weak for peak 
heights of the order of unity, never cancels out. Over- 
all, retaining the density correlation £(r) solely is not a 
reasonable approximation when the peak height does not 
exceed v <4. Although the exact value of the bias pa- 
rameters changes somewhat with the smoothing scale R / , 
their global behaviour varies little as 7 weakly depends 
on the filtering scale. Therefore, the above statements 
hold regardless of the exact amount of smoothing. 



F. Peak biasing : nonlinear and local ? 



Higher derivatives of the density field may also contribute 
to this general expression. However, (nonlocal) integrals 
of the linear density correlation are expected only in the 
evolved matter distribution when the non-Gaussianity in- 
duced by gravitational clustering is significant. 

Finally, it is worth noticing that, upon Fourier trans- 
formation, the peak power spectrum reads 



Pvk{v,k) = 6p k 



a\bl 



k 2 + 



P(k), (41) 



where 6 p k = b v /<jQ and P(k) is the power spectrum of 
the smoothed density field. The exact amount of scale- 
dependence induced by the nonlinear bias depends upon 
the exact value of v and R / . We defer a thorough inves- 
tigation of this effect to a future work. 



IV. CLUSTERING OF DENSITY PEAKS IN 
GAUSSIAN INITIAL CONDITIONS 

After a brief discussion on the peak-background split, 
we focus on the acoustic signature in the 2-point corre- 
lation of density maxima. We examine how the baryon 
acoustic oscillation changes with the filtering, the thresh- 
old height v and the small-scale behaviour of the trans- 
fer function. We find that the extra contributions 6,S(r) 
and b 2 ip(r) to the linear relation £ p k(^, r) = b 2 ^(r) can 
boost significantly the contrast of the acoustic peak. 



Equation (j3"6")l clearly differs from the linear, local re- 
lation £ p k(^, r) — b 2 ^(r) that would be expected if the 
peak overdensity Srip^x) = n p k(a:)/(?Jpk) were related 
to the underlying density field through the linear map- 
ping Sn p k(x) = b v v{x). However, it is worth noticing 
that eq. f|36|) is compatible with a nonlinear, local, deter- 
ministic biasing relation involving a differential operator. 
Namely, it can be explicitly checked that 

Sn pk (x) = b v v(x) + bQu(x) , (39) 

where u(x) = —V 2 8(x)/(j2 : leads to the correlation func- 
tion (|36|) . This demonstrates that, at large distances, 
£ p k(z^, r) can be thought as arising from a specific case of 
nonlinear local bias. We will see later (Sec. [Vj that this 
local mapping is also consistent with the peak pairwise 
velocity at first order. 

To make connection with the formalism introduced 
by [25j], we may conceive of a Taylor series 

i=0 ' i=0 

to describe the properties of the peak distribution at all 
separations and filtering scales. An expansion of the peak 
correlation £ p k(*/, f) beyond leading order will be required 
to determine the values of the b$ and b^f 1 when i > 1. 



A. Filtering scale and peak height 

The peak height v and the filtering radius Rf could in 
principle be treated as two independent variables. How- 
ever, in order to make as much connection with dark 
matter haloes (and, to a lesser extent, galaxies) as pos- 
sible, we will follow the Press-Schechter prescription [5lJ 
which is based on the critical density criterion issued from 
the spherical collapse dynamics [52| . Namely, we assume 
that density maxima with peak height v = 8 sc (z)/ao(Rf) 
identified in the primeval density field smoothed at scale 
Rf are related to dark matter haloes of mass Mf collaps- 
ing at redshift z. Moreover, we will only present results 
at redshift z = 0, at which the linear critical density for 
(spherical) collapse is 5 SC — 1.673, and the characteristic 
mass for clustering is M+ ps 3.5 x 10 12 M Q /h. While there 
is a direct correspondence between the massive cluster- 
sized haloes in the evolved density field and the largest 
maxima of the initial density field, it is unclear the ex- 
tent to which galaxy-sized haloes trace the initial den- 
sity maxima [7l| . For this reason, we will only consider 
mass scales in the range Mf ;>M*(0) or, equivalently, a 
smoothing radius Rf ~ 1 /i _1 Mpc. 

We note that the spherical infall model provides a local 
approximation to the collapse of a perturbation. How- 
ever, in the Press-Schechter approach, it is applied to 
random points in space and leads to linear local biasing 
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at large scales [53|, |5J] while, in the present work, it is 
applied to density maxima and leads to the specific type 
of nonlinear local biasing exemplified by eq. (|39|) . 

B. Peak-background split and the halo multiplicity 
function 

Before illustrating the impact of derivatives of the den- 
sity field on the baryon acoustic signature, we note that, 
in the limit v 3> 1, the peak correlation is amplified by 
an effective bias b 2 which is significantly smaller than 
the value v 2 /S sc derived for thresholded regions [3(j. As 
recognised in BBKS, this difference arises from the corre- 
lation between the peak height v and the peak curvature 
u. More precisely, in the spherical infall model, the linear 
Lagrangian bias 6p k = £ p k/(o"o£) of high density peaks 
that are collapsing at redshift z evaluates to 

fopk « -jj— (42) 

Osc 

in the limit r 3> 1. This should be compared to the 
expression derived in [55[ from the Press-Schechter 
formalism HU , 

v 2 - 1 

Omw = — ; , (43) 

o sc 

In this second approach, the clustering of haloes is de- 
scribed by the properties of regions above a given density 
threshold. In both cases however, the Kaiser limit v 2 /5 SC 
is recovered. This, however, does not apply to the bias 
factor derived by [5J| using the ellipsoidal collapse, 

av 2 -l /AAS 
&ST « —. , (44) 

Osc 

where a ~ 0.7. Assuming the peak-background split 
holds (30j . these various bias parameters predict mul- 
tiplicity functions vf(v) that have quite a different 
behaviour in the limit of large threshold heights. In par- 
ticular, the Sheth-Tormen (ST) multiplicity function is 
proportional to v exp(— av 2 /2) [561 ] . and exponentially 
deviates from the scaling inferred from o p k and Omw> 
which is vf(v) cx is 3 exp(— v 2 /2) and c>c i/exp(— v 2 /2), 
respectively. It is worth emphasising that the factor 
a = 0.707 was essentially determined by the number of 
massive haloes in the GIF simulations [57| and, therefore, 
is not a direct outcome of the ellipsoidal collapse dynam- 
ics. In fact, there is no compelling theoretical reason for 
a halo mass function whose high-mass end deviates ex- 
ponentially from the scaling exp(— v 2 /2). Furthermore, 
recent lines of evidence suggest that the high-mass tail, 
while being above the Press-Schechter (PS) mass func- 
tion [5l|, may depart from the Sheth-Tormen scaling 

In our opinion, it is likely that the true multiplicity 
function scales as exp(— v 2 /2) in the limit of large v. This 




50 100 150 

r / h _1 Mpc 



FIG. 4: The peak correlation £ p k(^, r) (solid curves) for 
three different smoothing lengths Rf = 2, 4 and 6 /i _1 Mpc 
(from bottom to top). These correspond to a mass scale 
Mf = 9.5 x 10 12 , 7.6 x 10 13 and 2.6 x 10 14 M Q /h, respec- 
tively. A peak height v — 5 sc /cro is adopted and yields the 
values v = 1.40, 2.15 and 2.88, respectively. The density cor- 
relation <7of(r) is plotted as the dotted-dashed curve. Dash- 
ing indicates negative values. The acoustic signature in the 
peak correlation depends on the threshold height v through 
the bias parameters b v , b v and 6^. Results are shown for the 
ACDM cosmology. 

would lead to a different parametrisation of the halo bias 
and mass function. Given the lack of a convincing physi- 
cal description of these quantities, one may, for instance, 
consider a phenomenological bias of the form 

fo L = ^ U -ci + -^—] (45) 

which, for a peak-background split, leads to a multiplicity 
function 

^)oc(l + ^)° 2/2PC3 ^e-^. (46) 

For ci <^3 and C3 ~ Oijc\ (which guarantees &l ~ in 
the limit v — ► 0), the biasing (|45| closely follows the 
peak scaling eq. (|42[) at large mass and, simultaneously, 
exhibits an upturn at low mass. Unfortunately, such a 
bias cannot be derived from an excursion set approach 
(upon which PS and ST are based), where c\ = 1 in- 
variably. This issue, which lies beyond the scope of the 
present paper, will be examined in a separate paper. 

C. Baryon acoustic signature 

We now turn to the behaviour of the peak correlation 
function. ^ p k(^, r) is shown in Fig. 0] for a filtering length 
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Rf = 2, 4 and 6 ft _1 Mpc. The mass enclosed in the 
Gaussian window thus is Mt = 9.5 x 10 12 , 7.6 x 10 13 and 
2.6 x 10 14 M Q /h, respectively. To illustrate, we have 
adopted a peak height v = S sc (z = 0)/oq such that 
v = 1.4, 2.1 and 2.9, respectively. In the spherical in- 
fall dynamics, a top hat overdensity enclosing a similar 
amount of mass would collapse at redshift z ~ 0. Fur- 
thermore, the density correlation <jQ^{r) is also shown for 
comparison as the dotted curve. 

The three correlations considered here exhibit a very 
different behaviour that reflects the strong dependence 
of the bias factors b„, b v and on the threshold height 
(see Sec. HUE]) . In particular, we find b v = -0.057, 0.847 
and 1.771 with increasing smoothing radius. As a result, 
for Rf = 2 /i _1 Mpc, the contribution of the term b^ip(r) 
dominates the others and strongly suppresses the ampli- 
tude of Cpkiv, r) relative to that of the density correla- 
tion. This term has the sign of ip{r) and features several 
oscillations across the BAO scale (see Figures [1] and [2]) . 
However, for peaks of threshold height v = 2.1 identified 
at smoothing scale Rf — 4 /i _1 Mpc, b^-ip(r) merely con- 
tributes to decrease the level of the minimum at distance 
r ~ 90 — 95 /i _1 Mpc. Interestingly, the term tiLT<(r) 
boosts significantly the contrast of the acoustic peak. 
This effect is still present, albeit weaker, for v = 2.9 
and Rf = 6 /i _1 Mpc. We also note that zero-crossings 
of £ p k(^, r) do not generally coincide with those of £(r), 
in agreement with numerical studies of the clustering of 
density maxima [3(1 [37j • 

Fig. E] further illustrates the sharpening of the acoustic 
peak due to correlations among derivatives of the density 
field. The density and the peak correlations are com- 
pared in the neighbourhood of the acoustic feature for 
the smoothing radii Rf — 4 and 6 /i -1 Mpc considered 
above. To emphasise the contrast of the acoustic peak, 
all the correlations have been rescaled such that, at a 
distance r = 70 /i~ 1 Mpc, their amplitude is equal to 
3. Fig. [5] nicely demonstrates the large impact of 6 2 E, 
which fully restores the acoustic signature of £,p\t(v,r) 
otherwise smeared out by the large filtering. The con- 
trast of the acoustic peak can even be enhanced relative 
to that of the unsmoothed (Rf = 0.1 h~ 1 Mpc) linear 
density correlation (dotted-dashed line). The effect is 
strongest for the density peaks identified at the smaller 
smoothing, Rf = 4 /i _1 Mpc. For these maxima, the dif- 
ference between the height of the (negative) minimum at 
r ~ 90 /i~ 1 Mpc and the maximum at ~ 105 h~ 1 Mpc is 
twice as large as in the linear density correlation. The 
enhancement is somewhat smaller, roughly 20 per cent, 
for the peaks at the filtering scale Rf = 6 h~ 1 Mpc. 
This shows that density maxima behave rather differ- 
ently than linearly biased tracers of the density field, 
whose acoustic signature cannot be larger than that of 
the linear matter correlation [l3l |. 

We now concentrate on the vertical lines which indicate 
the position of the local maximum. On the one hand, the 
top panel shows that smoothing in the density correla- 
tion generates a shift towards smaller scales, because the 
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FIG. 5: A comparison between the density correlation f(r) 
(top panel) and the peak correlation £pk(zA r) (bottom panel) 
around the BAO. The density field is smoothed with a Gaus- 
sian filter of width Rf = 4 and 6 /i -1 Mpc. The correspond- 
ing value of peak height is v = 2.1 and 2.9, respectively. For 
clarity, all the correlations have been rescaled such that, at 
separation r — 70 h~ Mpc, their amplitude is equal to 3. 
Also shown as the dotted-dashed line is the (unsmoothed) 
linear matter correlation. The vertical dashed lines indicate 
the position of the local maximum. The presence of 6^E(r) in 
the peak correlation restores, and even amplifies the acoustic 
peak otherwise smeared out by the large filtering. 6^E(r) also 
acts to reduce the shift induced by the smoothing. Results 
are shown for the ACDM cosmology. 



acoustic feature is not quite symmetric around its maxi- 
mum [14] . On the other hand, the presence of 6 2 E in the 
peak correlation acts in the opposite sense and compen- 
sates for the shift induced by the smoothing. We find the 
maximum to be close to its linear value ~ 105.0 /i _1 Mpc 
in both cases. More precisely, there is a small shift of 
<0.4 per cent towards larger scales. 



D. Sensitivity to the filter shape and the transfer 
function 

As discussed in BBKS, the filtering of the density field 
is an essential operation for power spectra covering a wide 
range of wavenumbers. However, the optimal choice of 
filter is disputable. Furthermore, the (analytic) proper- 
ties of the filtered density field can depend significantly 
upon the amount of power in small-scale fluctuations. 
It is, therefore, important to assess the influence of the 
smoothing operation and the small-scale transfer func- 
tion on the baryon acoustic signature in the correlation 
of density peaks. 
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To this purpose, we have repeated the numerical cal- 
culation of £ p k(^, r) using a top hat filter. To avoid diver- 
gence of the spectral moments and the correlation func- 
tions, we have introduced a high-fc cutoff whose func- 
tional form is motivated by the damping of fluctua- 
tions due to the free streaming of the dark matter par- 
ticle^). So far, we have considered a CDM cosmology 
in which the velocity dispersion of the dark matter par- 
ticle is negligible. By contrast, in Warm Dark Matter 
(WDM) cosmologies, the dark matter candidate(s) can 
suppress the matter power spectrum on galaxy scales 
r ~ 0.1 h~ 1 Mpc [59[ . The latter can be approximated 
as Pwdm(^) = T 2 (k)PcDM(k), where the transfer func- 
tion that accounts for the free-streaming cutoff has the 
form 60] 



T(k) 



1 



(ak) 



J 2p 



-5/p 



(47) 



Here, p w 1.12 and a depends upon the properties of 
the dark matter particles. Typically, 0.01 <a <j0.1 for 
thermal relics of mass ~ 1 — 10 keV. 

In spite of its compactness, the top hat filter has some 
inconvenient. Firstly, because of its slowly decaying tail, 
it produces a density field that is not differentiable for 
generic CDM power spectra [$fj. Secondly, it is good at 
discriminating peaks from the background field so long 
as the height of the latter is small, namely, when the 
background field is uncorrelated over scales comparable 
to the filtering length. By contrast, the Gaussian filter 
is less sensitive to high frequencies and thus fares better 
at picking up smoother objects. Indeed, the "true" filter 
may lie between these two extremes [6l[ . Notice that the 
sharp fc-space window will not be considered here as it 
leads to undesirable oscillations at all separations. 

Fig. [5] shows the baryon acoustic peak in the correla- 
tion of density maxima for the smoothing radii used in 
Fig. d] and O Note, however, that the filter mass scale is 
now roughly four times smaller than with the Gaussian 
window. The peak correlation is plotted for two values of 
the free-streaming cutoff, a — 0.01 and 0.1 (top and bot- 
tom panels). Also shown in both panels for comparison 
is the linear matter correlation (dotted-dashed line). At 
smoothing length Rf = 4 and 6 /i _1 Mpc, the enhance- 
ment of the acoustic peak is very significant for a = 0.1 
while, for a = 0.01, it is only 10-20 per cent. The main 
reason is a sharper power spectrum, which leads to a 
larger contribution of the correlations S(r) and ip(r) to 
£pk(v,r). At Rf = 4 /i _1 Mpc for instance, the spectral 
width is 7 = 0.48 and 0.26 for a = 0.1 and 0.01, re- 
spectively. This difference mostly arises because of the 
second spectral moment, which increases from 02 = 0.40 
to 0.76 upon the decrease in the free-streaming scale. Yet 
another interesting feature of Fig. [5] is the rather broad 
acoustic peak at filtering scale Rf — 2 /i _1 Mpc (see bot- 
tom panel), for which v = 0.96. This broadening follows 
from the fact that b n is negative at that value of thresh- 
old height. As a consequence, the oscillatory pattern of 
S(r) across the BAO (see Figures [1] and [2]) smears out 



the acoustic feature in ^ p k(i / , r). As seen from Fig- [3] this 
damping always occurs at sufficiently low values of the 
threshold height, v <1 regardless of the filtering length. 
It should also be noted that, unlike the correlation of 
density maxima, the BAO in the smoothed linear matter 
correlation is weakly insensitive to the small-scale be- 
haviour of the power spectrum. In ^ p k(z / , f) however, the 
BAO acquires an extra dependence upon the high-/c tail 
of the transfer function through the correlation functions 
E(r) and ip( r )- 
To summarise, 

• Both S(r) and ip(r) contribute to the correlation 
of density maxima and can affect the shape of the 
baryon acoustic signature for peak heights v <4. 

• ip(r) makes a significant contribution only in the 
range 1 <v <2, where b v and b v are much less than 
unity. 

• The contribution of £(r) increases with the spectral 
width 7. At constant filtering length, it increases 
with the amount of power suppression due to the 
small scale free streaming. 

• b v is positive (negative) for v J>1 {y <Jl). As a re- 
sult, the baryon acoustic peak is generally enhanced 
in £pk(^, r) when v >1, and damped out when v <1. 

These results depend upon the exact shape of the fil- 
ter and the transfer function. Clearly however, the effect 
cannot be reduced to a simple rescaling of the linear mat- 
ter correlation. This is due to the peculiar type of nonlin- 
ear local biasing, eq. (|39p . which involves the Laplacian 
of the density field. 



PAIRWISE VELOCITY OF DENSITY 
MAXIMA 



Thus far, we have explored the BAO signature in the 
correlation of maxima of the primordial density field. 
However, pairwise motions caused by small- and large- 
scale structures, redshift space distortions etc. are likely 
to degrade the acoustic signature, leading to a broadening 
and, possibly, a shift of the acoustic peak. A thorough in- 
vestigation of these effects is postponed to a subsequent 
work. Here, we consider a simple model in which the 
peak centres evolve according to the Zeldovich ansatz. 
This allows us to calculate the peak pairwise velocity at 
leading order, which is the main result of this Section. 
We show that, at first order, the peak mean streaming is 
consistent with the nonlinear local bias found in Sec. IIVI 
Dynamical evolution is also briefly addressed using the 
pair conservation equation. 
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FIG. 6: The correlation of density maxima that trace the 
density field smoothed with a Top Hat filter. The smooth- 
ing radii Rf — 2, 4 and 6 /i _1 Mpc correspond to a mass 
Mf = 2.5 x 10 12 , 2 x 10 13 and 6.8 x 10 13 M /ft, respectively. 
Results are shown for a WDM power spectrum with a cutoff 
scale a = 0.01 and 0.1 /i _1 Mpc (see text). The peak height 
is chosen such that v = S sc /ao, as before. In both panels, 
the dotted-dashed curve is the linear matter correlation. The 
peak correlation £ p k(^, r) for a = 0.01 and Rf — 2 is not 
shown as it is too much affected by numerical noise. 



FIG. 7: The correlations that contribute to the leading or- 
der mean streaming of peak pairs, equation (|51|l . These 
are compared to the line of sight pairwise velocity V(r) of 
ambient field points. Results are shown at a filtering scale 
Rf — 5 /i _1 Mpc. The correlation V — 7uE is strongly damped 
on scales less than the characteristic inter-peak distance tx Rf 
but, at large distances, it is unaffected by small-scale exclu- 
sion effects and closely follows the (scaled) mean streaming 
of random field points. The correlation <S — 7^11 can signifi- 
cantly contribute to small-scale streaming motions when the 
peak height is v <^3 (so that ~^b v ). 



A. Zeldovich approximation 

The Eulerian comoving position and proper velocity of 
a density peak can generally be expressed as a mapping 

x P k = q+ S(q,a) , v pk = aS(q, a), (48) 

where q is the initial position, S(q, a) is the displace- 
ment field and a is the scale factor. A dot denotes a 
time derivative. At first order, the pe ak p osition is de- 
scribed by the Zeldovich approximation [62| , in which the 
displacement factorises into a time and a spatial compo- 
nent, 

S = -L>(a)V$(g), S = -/3(a)V$(g) . (49) 

Here, $(q) is the perturbation potential linearly ex- 
trapolated to present time. Explicitly, $(q) = 
4>(q, a) /4:TrGp m (a)a 2 D(a) where <f>(q, a) is the Newtonian 
gravitational potential, p m (a) is the average matter den- 
sity and D{a) is the growth factor. (3{a) = HDf is pro- 
portional to the logarithmic derivative / = dlnD/dlna, 
which scales as f(a) ~ f2 m (a) ' 6 for a wide range of CDM 
cosmologies [63|]. Finally, H(a) is the Hubble constant. 

Such a simple model cannot account (among other 
things) for the internal properties of peaks 64]. Fur- 
thermore, it provides a very limited description of the 



late-time distribution of density maxima such as cluster- 
or galaxy-size haloes [23l. [65|. Notwithstanding this, it 
is not intended to be realistic, but only to capture the 
weakly nonlinear regime reasonably well. A more sophis- 
ticated approach can be found in [66J for instance. 

The peak pairwise velocity, or mean streaming (67I . |68| , 
is now obtained from the statistics of the (proper) matter 
velocity field v = — af3(a)VQ(q). The complication arises 
from the fact that the latter has to be evaluated at those 
maxima of the density field. 



B. Mean streaming of peak pairs 

We introduce the normalised velocity field u(q) = 
v(q) / (af3(a)<T-i) for subsequent use, and define uia(r) as 
the average number weighted pairwise velocity [u(a^) — 
u(cci)] • r along the line of sight. 

The calculation of the peak pairwise velocity is more 
intricate than the peak correlation since we have three 
additional degrees of freedom. Nevertheless, it closely 
follows the analysis described in Sec. Mil Details of the 
calculation can be seen at Appendix |XJ The peak mean 
streaming weighted by the number density of peaks at q 1 



14 



and q 2 eventually reads 

(l+fpk)ui2 = [M^.7) + M^,7)](V-7«S) (50) 
+ [6 c (^i,7) + 6 c (iv 2 ,7)](<S-7„n) . 

For comparison, the pairwise velocity of the matter dis- 
tribution is [1 + crQ^(r)]ui2(r) — 2<r V(r). The first term 
in the right-hand side of eq. (f5Tj) is similar to the mean 
streaming derived for locally biased tracers of the density 
field [69(. The second term arises because of the partic- 
ular nature of the bias of density maxima. Interestingly, 
eq. (|5Tj) can again be thought of as arising from the non- 
linear local bias eq. (|3"9")) if we choose the peak velocity 
field to be 

u pk (a;) = u(x) - -f v ri{x) . (51) 

In this continuous approach, the peak velocity field is still 
unbiased with respect to the matter velocity field u(x), 
but it receives a contribution from the first derivative of 
the density, r)(x) = V8(x)/a\, that is proportional to ~/ v . 

The sign and the strength of the peak pairs flow depend 
upon the detailed behaviour of the functions V — 7„3 
and S — j v Tl. As seen in Fig. [7] where Rf = 5 /j _1 Mpc 
for illustration, the former is negative at all separations 
r <200 /i _1 Mpc. By contrast, the latter is positive at dis- 
tances larger than a few smoothing radii but goes neg- 
ative at smaller scales, regardless of the exact value of 
Rf. The mean streaming V(r) of random field points 
is also shown in Fig. [7J for comparison. Peak-peak ex- 
clusion leads to a deficit of pairs at separation r <R* 
comparable to the filtering scale. It adds to the smooth- 
ing and further damps the relative velocity of peaks out 
to distances that are much larger than the typical ex- 
tent ~ Rf of density maxima. This is the reason why 
the correlation V — 7^S is strongly suppressed relative to 
V(r) when r <30 /i~ 1 Mpc. Still, the term proportional 
to S — j v Il is most strongly negative at distances of the 
order of the smoothing length and, therefore, could re- 
store significantly the small-scale mean streaming when 
the threshold height is less than v <3 (for which >b v ). 
At large enough separations r >50 h~ 1 Mpc however, the 
mean streaming of peak pairs is unaffected by small-scale 
exclusion effects and closely tracks the pairwise velocity 
V(r) of ambient field points. 

To illustrate the impact of the correlation function 
S — 7^11 on the mean streaming, we show in Fig. [5] the 
peak pairwise velocity for the smoothing radii Rf = 2, 4 
and 6 h~ 1 Mpc considered in Sec. lIVI (recall that the peak 
height is specified by the relation v = 1.673/ ao(Rf)). 
Also shown for comparison is the mean streaming of 
the (unsmoothed) matter density field (dotted-dashed 
curve). At the smallest filtering scale for which v = 1.4, 
the bias parameters are — b v » 0.06 <C w 2.32 so 
that S — 7„II is the dominant contribution at all sepa- 
rations. The resulting strong "inward" transport at dis- 
tances less than a few /i _1 Mpc reflects the fact that 
these small peaks tend to accrete onto high density re- 
gions. At separation r ~ 10 h~ 1 Mpc, there is a positive 
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r / h"'Mpc 

FIG. 8: The mean streaming of peak pairs for density maxima 
identified at smoothing scale Rf = 2, 4 and 6 ft. - Mpc and 
with peak height v — 1.40, 2.15 and 2.88, respectively (as in 
Fig. [4]| . The dotted-dashed curve shows the pairwise velocity 
of the (unsmoothed) underlying density field. The smallest 
peaks tend to accrete onto high density maxima. However, 
they move apart from each other relative to the matter dis- 
tribution due to peak-peak exclusion. 

net flow presumably owing to the fact that the peaks 
fall onto nearby distinct overdense regions. Notice that, 
at all separation, the mean streaming of these maxima 
is larger than that of the density field, indicating that 
these small peaks move apart from each other (in an av- 
erage sense) relative to the matter distribution. At larger 
filtering scales, the contribution of the first term in the 
right-hand side of eq. (f5Tj) increases with the smoothing 
length as seen from the progressive disappearance of the 
broad bump at r ~ 10 /i _1 Mpc. The maxima identified 
at Rf = 6 /i _1 Mpc [y = 2.9) stream towards each other 
relative to the underlying density field, but their relative 
motion is strongly suppressed at distances r <10 /i -1 Mpc 
due to the exclusion effect mentioned above. 



C. Pair conservation equation 

So long as peaks do not merge, the time evolution of 
the peak correlation function £ p k is governed by the pair 
conservation equation. Transforming the time variable 
to the scale factor, this equation can be written as 

^ = -Df<r-i~\r 2 (l + ^)u 12 (r)] , (52) 

where r and u\2(r) are the comoving separation and 
scaled pairwise velocity, respectively. The root-mean- 
square (rms) variance cr_i (computed from eq. [T|) defines 
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the length scale <r_i w 9.2 /i _1 Mpc for the values of cos- 
mological parameters used here. 

Following the approach outlined in pjj, the general 
solution of eq. (|52|) can be found by solving the charac- 
teristic equation 

^ = Df<j- lUl2 (r) . (53) 
da 

This equation gives r(a) which, upon insertion into the 
pair conservation equation ([52"]) , allows us to write down 
a first order ordinary differential equation along the char- 
acteristics, 

dln[l + g pk (r, a )] nf 1 d[r 2 u 12 (r)) 
da r z or 

where it is understood that r = r(a) [l7j . 

We will not attempt to solve eq. [54] since, as recog- 
nised by [lil nonlinearity in the divergence of the 
pairwise velocity, which is lacking here, is a crucial in- 
gredient in the redshift evolution of the baryon acoustic 
signature. Instead, we will simply estimate the first or- 
der change in the initial separation of peak pairs, Ar , 
induced by coherent motions across the acoustic scale, 
ro ~ 105 ft. Mpc. To proceed, we assume that the peaks 
move according to the Zeldovich ansatz described above, 
an approximation expected to be valid only in the early 
(quasi-linear) stages of gravitational clustering. Owing to 
the near constancy of the peak pairwise velocity at those 
scales, we can write Aro » a^iu^iro) / Dfda where 
J Dfda as 0.56. For the maxima considered above, we 
find Ar (r = 105) = +0.010, -0.25 and -0.74 /i^Mpc 
with increasing Rf. For comparison, A(tq = 105) w 
0.21 /i -1 Mpc for the dark matter. These values are con- 
sistent with those found by [ljj. Therefore, at the linear 
order, changes in the acoustic signature of £ p k(?") will be 
roughly at the percent level. This suggests that some 
of the enhancement of the BAO in the initial correla- 
tion of density maxima may survive in the correlation of 
high redshift density peaks. Clearly, a thorough numeri- 
cal investigation and detailed analytic modelling will be 
needed to ascertain how much of this effect propagates 
into the late-time clustering of galaxies. 

VI. CONCLUSIONS 

We have investigated the strength of the baryon acous- 
tic signature in the 2-point correlation of maxima of the 
linear (Gaussian) density field 5(x). To this purpose, we 
examined in Sec. IIIII the large-scale asymptotics of the 
peak correlation £ p k(»") and derived the leading order con- 
tribution, eq. (f3T>)) . In contrast to the analysis of BBKS, 
spatial derivatives of the linear density correlation £(r) 
were included in our derivation. These derivatives are 
not negligible for generic CDM power spectra, especially 
around the BAO scale where they exhibit large oscilla- 
tions, we find that the leading asymptotic behaviour of 



the peak correlation is governed by three terms : a term 
previously derived in BBKS plus two terms involving the 
spatial derivatives S(r) and ip(r) of the linear density 
correlation. The relative contribution of these functions 
is controlled by two independent bias parameters, &„ and 
&f (The third being b v = 2"fb u bQ, see eq. [37]). We also 
showed that the large-scale asymptotics of £ p k(?") can be 
thought as arising from a nonlinear, local biasing rela- 
tion, eq. 1391 involving second derivatives of the density 
field. 

In Sec. IIV] we demonstrated that those extra terms 
have a large impact on the correlation of density maxima 
in the vicinity of the BAO. The results are sensitive to 
the exact value of the threshold height v, the smooth- 
ing length Rf, the filter shape and the high-fc tail of 
the transfer function. For the Gaussian filter adopted 
throughout this paper, the contrast of the baryon acous- 
tic signature can be significantly enhanced relative to 
that in the linear matter correlation when the peak height 
is in the range 1 <v <3. This boost originates from the 
oscillatory behaviour of S(r) and ip{r) around the sound 
horizon scale. For instance, we find that, at filtering scale 
Mf = 8x 10 13 Mts/h, the contrast of the BAO in the 
correlation of density maxima of height is about twice as 
large as in £(r). The amplification fades as we go to larger 
peak height. For a peak height of the order of unity, 
Cpk(^) can exhibit several bumps which reflect those of 
ip(r) around the BAO scale. For a threshold height less 
than <1, the original acoustic peak is smeared out by 
the negative contribution of the term 6 I? E(r). 

To avoid the divergence of the (fourth order) spatial 
derivative ip(r) of the density correlation, we have filtered 
the density field with a Gaussian window. The main 
drawback of this window function is the lack of a well- 
defined mass and spatial extent associated to the density 
fluctuations. A top hat filter appears better motivated in 
the context of, e.g., the spherical infall model, although it 
does not produce an infinitely differentiable density field 
for generic CDM power spectra. Furthermore, the differ- 
entiability of the density field depends strongly upon the 
small-scale behaviour of the transfer function. Fluctua- 
tions in the matter density are damped on scales smaller 
than the free-streaming length of the dark matter parti- 
cle. In the CDM cosmology considered here, the velocity 
dispersion of the dark matter particle is negligible. By 
contrast, Warm Dark Matter (WDM) particles such as 
massive neutrinos for example can suppress the matter 
power spectrum on galaxy scales 59]. For these reasons, 
we also discussed in Sec. IIVI how the BAO changes with 
the window function and the free-streaming cutoff. We 
found that the correlation of density maxima is more 
sensitive to the properties of the dark matter particle(s) 
than the matter correlation itself. This follows from the 
dependence of the peak biasing upon the second deriva- 
tives of the density field. However, whether the baryon 
acoustic signature in the clustering of peaks varies sig- 
nificantly with the nature of dark matter remains to be 
determined. 
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In Sec. El we calculated the pairwise velocity of peak 
pairs at leading order. We showed that it is consistent 
with the nonlinear local biasing relation inferred from 
the 2-point correlation of density maxima, provided that 
the peak velocity field receives a contribution from the 
gradient of the density field. Explicitly, the leading-order 
peak correlation and mean streaming can be derived from 
the nonlinear local biasing relation 

5n p k(x) — b u 6(x) — b(;V 2 S(x) 
Upk(x) = u(x) - 7„V<S(x) , (55) 
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where we have dropped some factors for clarity and u(x) 
is the linear matter velocity field. This particular bias 
relation may be helpful to translate the Lagrangian anal- 
ysis performed in this paper into quantitative predictions 
for the baryon oscillation in the low redshift distribution 
of galaxies, which is currently the primary observable 
proxy of the baryonic acoustic oscillations. Using the 
formalism introduced by (25| . one may conceive of so- 
phisticated extensions of the halo model [70J that would 
include derivatives of the density field, so as to ascertain 
how much of the amplification of the acoustic signature in 
the initial clustering of density maxima propagates into 
the late-time correlation of galaxies and clusters. Ex- 
tensions which, as a general criterion, reproduce the ob- 
served properties of the galaxy distribution, would pro- 
vide an interesting complement to current local biasing 
models. 

We emphasise that the calculations presented in this 
paper are performed in the initial conditions. As non- 
linearities progress, the late-time acoustic signature is 
smeared out by structure formation as reported by many 
authors using N-body simulations 0]. This might ex- 
plain why numerical investigations of the clustering of 
dark matter haloes have not shown thus far any evidence 
for an amplification of the BAO. Interestingly however, 
preliminary results from a large suite of N-body simu- 
lations hint at a an enhancement of the contrast of the 
baryonic signature in the clustering of low redshift dark 
matter haloes [72[. In this regard, it is also worth notic- 
ing that the clustering of the SDSS (Sloan Digital Sky 
Survey) LRG (Luminous Red Galaxies) sample [2j], for 
instance, shows a slightly sharper acoustic peak than ex- 
pected from linear theory and smearing due to nonlin- 
earities. However, one should remember that the data 
points are strongly correlated, so that a very high acous- 
tic peak is actually allowed by the current ACDM cos- 
mology. Future redshift surveys such as ADEPT, BOSS, 
CIP, DES, HETDEXjLSST, Pan-STARRS, PAU, Wig- 
gleZ or WFMOS /73], which will obtain redshifts for 
millions of galaxies, should achieve an exquisite precision 
on the shape of the baryon acoustic signature in the clus- 
tering of galaxies. Beyond the nature of dark energy, a 
precise measurement of the BAO could also place con- 
straints on galaxy biasing and the physical mechanisms 
that cause it. 



APPENDIX A: MEAN STREAMING OF PEAK 
PAIRS 

1. Correlations of velocity field 

Let us introduce the scaled velocity field u{q) = 
v(q)/(a/3(a)er_i). The rms variance a(3a-i is the three- 
dimensional proper velocity dispersion of random field 
points, which is ~ 430 kms^ 1 at present time in the cos- 
mology considered here. Also, the notational shorthand 
Am will designate the difference u{q 2 )~u(q 1 ). The auto- 
correlation of the velocity and its cross-correlations with 
the fields 7/i, v and Qj can be written as 

(vi(q 1 )v j (q 2 )) =U 1 (r)f i f j + %(r) 6^ (Al) 

M 9i)^(92)) = 7u fi(*0 hrj + 7t, 6(r) 
(v(qi)Vi(q 2 )) = V(r)f t 

= Si(r)fifif m + S 2 (r) {hSi m + fid im + f m 6u) . 

Here Vi(q) designates the components of u(q). Notice 
that U\\(r) = Ui(r) +U 2 {r) and U±(r) = U 2 {r) are the 
radial and transverse correlation functions of the velocity 
field [73| . For sake of completeness, the various angle 
average correlations are 
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for the Gaussian density field considered here. The func- 
tions £1 and £2 satisfy the relation £i(r) + ^^(r) = 
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£(r). Like the spectral width 7, the parameter 7^ = 
(7q/((7_i(7i) characterises the range over which the veloc- 
ity power spectrum oc k~ 2 A 2 (k) is large. It should be 
noted that the latter peaks on scale much larger than 
the density power spectrum. Also, the correlation V(r) 
is proportional to the mean streaming of ambient field 
points, 

<[1 + 5{ qi )] [1 + S(q 2 )] Au ■ r) = 2a a V(r) , (A3) 
which is mass weighted by the densities at q 1 and q 2 . 

2. Mean streaming at leading order 

The calculation of the peak pairwise velocity is more 
intricate than the peak correlation since we have three 
additional degrees of freedom, but it closely follows the 
analysis described in Sec. IIIIl 

The line of sight pairwise velocity weighted over all 
pairs with comoving separation r can be expressed as 

[l + £pk(r)]uu(r) = (n pk y 2 (A4) 
x ^: /dO f .dy 1 dy 2 (An- r) n^q^n^^q^) P(y 1 , y 2 ; r) 

The local peak density n p k(g) is given by equation ((7|, 
supplemented by the appropriate conditions to select 
those maxima with a certain threshold height. To ob- 
tain the average pair velocity as a function of separation 
r, we need to calculate the 2-point probability distribu- 
tion for the variables y T — (^,77^,^,^)- At zero lag, 
both Vi and rji are uncorrelated with the density and the 
Hessian £4- Hence, the covariance Mi of the components 
(vi ,r]i,v) is a 7 x 7 block matrix which reads 



1/31 7 „/3I 3 xi 
Mi = I 7„/3I 1/31 3x i 

0lx3 0ix3 1 



(A5) 



Similarly, the covariance M2 of the Hessian, and the 
cross-covariance M3 between C,a and the entries (i>i ,^,1;) 
are 
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03x3 3X 3 -7/313x1 
03x3 03x3 03x1 



(A6) 



Proceeding as in Sec. IIIIl we now consider the regime 
where all the correlations are much less than unity. The 
2-point probability distribution P(y 1 , y 2 ; r) can thus be 
expanded in the small perturbation B(r), 

P{yx,V2\r) « (l + ^BM- 1 ^) P(y l )P{y 2 ) ■ 

(A7) 

Here P(y) designates the 1-point probability density. As 
before, the (now 13 x 13) matrix B(r) denotes the covari- 
ances at different comoving positions. It has a (unique) 



harmonic decomposition in term of the matrices B^" 1 
(equation (fT2"l) . The computation of these matrices is, 
however, unnecessary as we will see later. Furthermore, 
the quadratic form Q(yi,y 2 ) now reads 
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We note that the velocity dispersion of density maxima 
is lower by a factor 1 — 7^ than that of random field 
points [29]. One has j v « 0.43 for a smoothing length 
Rf = 5 /i _1 Mpc. Moreover, eq. (|A8[) leads to a one- 
point probability distribution P(y|peak) oc exp[— Q(y)] 



separable into the product Pu(ui)-P^(^, 00, where Pi 



is the one-point distribution of the density and its sec- 
ond derivatives, and P v (vi) is the velocity distribution of 
peaks, 

33/2 
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3t/ 
'2(1-72) 



(A9) 

The separability of the one-point distribution separabil- 
ity considerably simplifies the calculation. 

Taking the product (An • r) B(r) mixes the various 
multipole matrices B^' m , so that the result depends on 
the correlation functions of Ui, 77^, v and C,a m a rather 
complicated way. Averaging over the directions gives 

1 



B = — dn~ r (Au - r) M^BM" 1 

47T 



B3 B 2 

(A10) 

where the block matrices B^ have the same dimensions 
as Mi . The minus sign in the right-hand side of eq. IA10I 
arises from the negative parity of the correlations (vi Q m ) 
and (r)i(; m ). Owing to the angular average, the calcu- 
lation of the B^'" 1 can be avoided by writing down the 
entries of B^(r) using the relations eq. ([3]) and (|A3|) . and 
retaining only those components involving odd products 
of the unit vector ?v A tedious calculation shows that 
Bi(r) and B 3 (r) can be cast into the form 
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L>i = I 3x3 3x3 ~a 2 Au 
aiAvJ a 2 Au T 
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The functions a^(r) are 
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where we have omitted the explicit r-dependence of the 
correlations for brevity. The 3x3 matrices have the 
components Avi of the vector Ait as entries, 
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, S(r) = 


Si + 
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The matrix B 2 is identically zero. 

The rest of the calculation is easily accomplished owing 
to the factorisation of the one-point probability distribu- 
tion P(y|peak). Notice that the scalar yjBy 2 contains 
terms linear and quadratic in U\ and Uq. After integrat- 
ing out the velocities, the linear terms vanish and we 
eventually find 



d 3 tiid 3 t*2 yjBy 2 P(Mi|peak)F(tt 2 |peak) = (A15) 
ai (ui + v 2 ) + (7«! - a 3 ) (trCi + trCa)] (l ~ 7«) • 



Transforming to the set of variables (ui,Vi,Wi) and sub- 
stituting the expressions (|3T[) of the bias parameters b v 
and 6f, the mean streaming of peak pairs can be recast 
into the form of equation (l5Tj) . 
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